Prediction of nuclear reaction rates for astrophysics 
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The investigation of explosive nuclear burning in astrophysical environments is a chal- 
I lenge for both theoretical and experimental nuclear physicists. Highly unstable nuclei are 

produced in such processes which again can be targets for subsequent reactions. The major- 
Psj . ity of reactions can be described in the framework of the statistical model (compound nucleus 

mechanism, Hauser-Feshbach approach), provided that the level density of the compound 
nucleus is sufficiently large in the contributing energy window [||]. Among the nuclear proper- 
ties needed in this treatment are masses, optical potentials, level densities, resonance energies 
CS| ! and widths of the GDR. All these necessary ingredients have to be provided in as reliable a 

way as possible, also for nuclei where no such information is available experimentally. 
^ ■ A recent experiment |§] has underlined that the low-energy extrapolation of the widely 

CD ! used optical a+nucleus potentials may still have to be improved. Currently, there are only 

I few global parametrizations for optical a+nucleus potentials at astrophysical energies. Most 

■ global potentials are of the Saxon-Woods form, parametrized at energies above about 70 

, MeV, e.g. Q. The high Coloumb barrier makes a direct experimental approach very difficult 

at low energies. More recently, there were attempts to extend those parametrizations to 
energies below 70 MeV Q. Early astrophysical statistical model calculations ^ made use 
of simplified equivalent square well potentials and the black nucleus approximation. Improved 
calculations Q employed a phenomenological Woods-Saxon potential Q, based on extensive 
data However, it was not clear how well all these potentials would work for heavy targets 
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■ with ^4 > 60 or in the thermonuclear energy range. 
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Most recent experimental investigations [10, 11 1 found a systematic mass- and energy- 
dependence of the optical potentials and were very successful in describing experimental 
scattering data, as well as bound and quasi-bound states and B(E2) values, with folding 
I potentials. Based on that work, a global parametrization of the volume integrals can be 

found [^]. In this description, the real part of the nuclear potential is given by a folding 
potential Vf{r,E). The imaginary part W{r,E) is of Woods-Saxon shape with a strongly 
energy-dependent depth. Nuclear structure and deformation information determines the 



shape of the energy-dependence by including level density dependent terms [12|. 

It is easy to show that the final transmission coefficients are not only sensitive to the 
strength of the potential but also to its geometry. Experimental data seemed to indicate 
that the geometry may also be energy-dependent |l^. At low energies, the diffuseness 
of the standard volume Woods-Saxon potential had to be set to smaller values, while the 
radius parameter was increased, in order to be able to describe experimental scattering data. 
This can be understood in terms of the semi-classical theory of elastic scattering [0] which 
shows that the relative importance of contributions from different radial parts of the potential 
depends on the energy. It was shown |15| that the predicted enhanced surface absorption at 



low energies can be described by an increased surface Woods-Saxon term. Thus, the artificial 
change in geometry in the description of scattering data results from the use of a volume term 
only. Consequently, the optical potential proposed here contains an imaginary part which is 
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given by the sum of a volume Woods-Saxon term and a surface term: 

V{r,E) = Vc{r) + Vf{r,E) +i(w,{E)f{r,R,a) - Ws{E)-^f{r,R,a)) (1) 



with 



f{r,R,a) = \l + e^] \ W^{E) = C - ae'^^ , Ws{E) = D + -fe-^'^ . (2) 



The depths of the potentials are exponentially dependent on the energy, with the volume 
depth Wy increasing and the surface depth Wg decreasing when going to higher energies. 
An increasingly dominant surface term at low energies leads to similar effects as reducing 
the diffuseness of a pure volume Woods-Saxon potential. The coefficients are related to 
the height of the Coulomb barrier and the microscopic and deformation corrections as used 



in Ref. [^]. The total volume integral is still given by the relation derived in Ref. |12|. 
The energy-dependence of the ^'^^Sm(a,7)^^^Gd excitation curve Q at low energies can be 
reproduced by such a description. Nevertheless, more experimental data is needed which 
should be consistently analyzed at different energies with optical potential parametrizations 
similar to the one used in Ref. ||l5[| . 

Based on the well-known code SMOKER P], an improved code for the prediction of 



astrophysical cross sections and reaction rates in the statistical model has been developed |16|. 
Among other changes, it includes an improved level density description |||], updated data sets 
of experimental level information, as well as the new a+nucleus potential. It also allows to 
treat isospin effects which are especially important in a capture reactions on self-conjugate 
target nuclei and in proton capture reactions above the neutron separation energy. For a 



more detailed presentation of the code and possible isospin effects, see |16 |. 
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